;***************************************************************
;discription：计算植被层对入射光合有效辐射的吸收比例
;
;param1: nrows(input) 行
;param2: ncols(input) 列
;param3: ndvi(input)归一化植被指数  格式为矩阵
;param4: vegetation(input,option) 植被类型图
;param5: pfar(output) 入射光合有效辐射的吸收比例
;
;author: wangqiang( wensten )  2015-09-14  chengdu
;***************************************************************
PRO PRO_APAR_FPAR,nrows, ncols, ndvi, vegetaion=vegetation, pfar=pfar
  ;不同植被类型比值植被指数最大最小值可能不同;
  
  ;本程序采用 "中国陆地植被净初级生产力遥感估算"-朱文全等(2006,植物生态学报)论文中估算
  sr_min=1.05
  sr_max=4.46
  ndvi_min=0.023
  ndvi_max=0.634
  alpha = 0.5
  
  ;取固定值，与植被类型无关
  pfar_min=0.001
  pfar_max=0.95
  
  FOR yindex=0,nrows-1 DO BEGIN
    FOR xindex=0,ncols-1 DO BEGIN     
      ;获取ndvi像元值
      ndvi_cellvalue = ndvi[yindex,xindex]
      
      IF ndvi_cellvalue EQ -9999.0 THEN BEGIN
        ;无数据区
        pfar[yindex,xindex] = -9999.0
      ENDIF ELSE BEGIN
        ;归一化植被指数计算所得
        pfar_ndvi=((ndvi_cellvalue - ndvi_min)/(ndvi_max - ndvi_min))*(pfar_max - pfar_min ) + pfar_min
        
        ;比值植被指数计算所得
        sr = (1.0 + ndvi_cellvalue)/( 1.0 - ndvi_cellvalue)
        ;光合有效辐射的吸收比例
        pfar_sr = ((sr - sr_min)/(sr_max - sr_min))*(pfar_max - pfar_min ) + pfar_min
        
        ;加权平均
        pfar[yindex,xindex]= alpha*pfar_ndvi + (1-alpha)*pfar_sr
        
      ENDELSE
    ENDFOR
  ENDFOR
  
END

;
FUNCTION FUNC_FPAR_SRMIN,b1,b2
  ;get min image
  b_min=b1>b2
  return, b_min
END